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O ■ ABSTRACT 



Using stacked Sloan Digital Sky Survey spectra, we present the detection of 
[O II] AA3727,3730 nebular emission from the galaxies hosting Ca II absorption line sys- 
tems and galaxies hosting Mg ll-selected damped Lyman-a (DLA) absorbers. Both samples 
^ ■ of absorbers, 345 Ca II systems and 3461 Mg ll-selected DLA systems, span the redshift 

' interval 0.4 ^ Zabs < 1-3; all of the former and half the latter sample are expected to be bona- 

00 . fide DLAs. The measured star formation rate (SFR) per absorber from light falling within the 

' SDSS fibre apertures (corresponding to physical radii of 6 — 9/i^^ kpc) is 0.11 — 0.14 Moyr^^ 

for the Mg ll-selected DLAs and 0.11 - 0.48 Moyr"^ for the Ca ll-absorbers. These results 
represent the first estimates of the average SFR in an absorption-selected galaxy population 
^5 ■ from the direct detection of nebular emission. Adopting the currently favoured model in which 

DLAs are large, with radii > 9h^^ kpc, and assuming no attenuation of the [O ll] emission by 
dust, leads us to conclude that the SFR per unit area of Mg ll-selected DLAs falls an order of 
CIh' magnitude below the predictions of the Schmidt law, which relates the SFR to the H I column 

density at z^O. While DLA sightlines are known to contain little dust, the unknown geometry 
^ ■ of the dust distribution in the galaxies causes the main uncertainty in our results. 

' The contribution of both DLA and Ca II absorbers to the total observed star formation 

, rate density, p*, in the redshift range 0.4 < z < 1.3, is small, < 10% and < 3% respectively. 

The result contrasts with the conclusions of Hopkins et al. that DLA absorbers can account for 
the majority of the total observed p* in the same redshift range. The disagreement is a direct 
, consequence of the much lower SFR per unit area we observe than predicted by the Schmidt 

\^ ' law. Our results effectively rule out a picture in which DLA absorbers are the sites in which a 

\ large fraction of the total p* at redshifts z < 1 occurs. 
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1 INTRODUCTION 

Uncovering the physical properties of the host galaxies of strong 
quasar absorption line systems is a key topic in extragalactic astron- 
omy today. These absorbers are expected to select galaxies based 
on gas cross-section, rather than luminosity, providing a comple- 
mentary view of the galaxy population to that obtained from tra- 
ditional stellar luminosity selected galaxy samples. However, with 
the possible exception of the common Mg II AA2796,2803 absorp- 
tion line systems, the identification of large samples of absorption- 
selected objects with known emission properties has not been pos- 
sible, despite considerable effort. A complete picture relating the 
gas cross-section, ionisation state and chemical abundances of the 
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interstellar media and gaseous halos surrounding galaxies to the 
emission properties of the associated galaxy remains elusive. 

At high redshift {z>2), observational efforts have focused 
on the highest neutral hydrogen (H I) column density systems. 
Damped Lyman-a systems (DLAs). These systems contain the ma- 
jority of neutral gas at high redshifts and may provide the fuel 
reservoirs for subsequent star formation in galaxies (Lanzetta et al. 
1995). Emission, in general from Lyman-a, from less than a dozen 
host galaxies of high redshift DLAs has been detected and spec- 
troscopically confirmed: a summary is given in Weatherley et al. 
(2005). Adopting only modest corrections to allow for the presence 
of dust, both detections and limits to the inferred star-formation 
rates (SFR) of these systems are relatively large, >lOM0yr~^. 

The small, inhomogeneous samples of counterpart galaxies 
have precluded direct estimates of the total SFR density (p*) of 
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galaxies associated with liigh redsliift absorption systems; liowever, 
two methods have been used to obtain indirect estimates. Hopkins 
et al. (2005) extrapolated the low redshift relation between H i col- 
umn density and SFR to high redshift; while Wolfe et al. (2003a, 
2003b) focused on the use of absorption from the fine-structure lev- 
els of the ground state of C II to estimate the SFR per unit area 
in DLAs and to calculate their contribution to p* at high redshift. 
Both approaches involve significant assumptions and the conclu- 
sions of the two studies are only marginally consistent, with Wolfe 
et al. (2003) advocating a larger contribution to the SFR density 
than Hopkins et al. (2005). 

At lower redshifts {z<l), galaxies associated with absorption- 
line systems are more readily observed because the angular sepa- 
ration between the galaxy and quasar is greater and deep imaging 
allows the detection of intrinsically fainter galaxies. Deep imag- 
ing and follow-up spectroscopic surveys of the fields of both DLAs 
and Mg ii absorbers have revealed host galaxies with a wide range 
of morphological types and at varying impact parameters from the 
background quasar (Bergeron & Boisse 1991; Le Brun et al. 1993; 
Steidel et al. 1994, 1997; Le Brun et al. 1997; Rao et al. 2003; 
Chen et al. 2005). Small samples, with incomplete spectroscopic 
follow-up, again hinder inference of the SFR density of absorption- 
selected galaxies. For a recent review on the nature of galaxies as- 
sociated with Mg II absorbers see Churchill ct al. (2005). 

Recently, Wild & Hewett (2005a) identified a new class 
of quasar absorption line system selected by their strong Ca ii 
AA3935, 3970 absorption doublet^ from quasar spectra in the Sloan 
Digital Sky Survey (SDSS). These rare (number densities ~ 30% 
of DLAs) absorption line systems contain significant quantities of 
dust - {E(B—V))^Q.l - in contrast to other samples of quasar ab- 
sorption line systems (Pettini et al. 1997; Murphy & Liske 2004; 
York et al. 2006). The dust and metal line properties of the Ca II 
absorbers strongly suggests they have H I column densities greater 
than the nominal limit for DLAs of 10^°'^ atoms cm~^ (Wild et al. 
2006, hereafter WHP06). 

The presence of strong Ca ii absorption in absorption line sys- 
tems is of particular interest for studying the properties of the inter- 
stellar medium (ISM) of galaxies. Due to the convenient position 
of its resonance lines in the optical spectrum, Ca II has been well 
studied in the disk of the Milky Way (MW): its distribution is found 
to vary smoothly within the warm neutral medium and extend to a 
scale height of about 1 kpc from the plane of the disk (Edgar & 
Savage 1989; Welty et al. 1996; Welsh et al. 1997; Hunter et al. 
2006). Calcium is severely depleted onto dust grains in the MW 
(Savage & Sembach 1996) and, due to its second ionisation poten- 
tial of 11.87eV being below the ionisation potential of hydrogen, 
Ca II is a trace ionisation state in the ISM of galaxies - most cal- 
cium is in the form of Ca III. The smoothness of the distribution, 
despite the sensitivity of gas phase Ca ii to ionising radiation fields 
and dust, have led to the suggestion that it traces warm neutral gas 
clouds within the MW, along with similarly depleted elements such 
as titanium (Hunter et al. 2006). Alternative explanations for the 
detection of strong Ca ii lines in absorbers at z^l could be the de- 
struction of dust grains, perhaps by shocks releasing calcium into 
the gas-phase, or lower radiation fields, perhaps within high volume 
density clouds which afford a degree of self shielding. 

WHP06 investigated the dust content of the absorbers by mea- 
suring both the reddening of the background quasar and dust de- 
pletion patterns of different ions detected in stacked SDSS spectra. 



Vacuum wavelengths are used throughout this paper. 



A dust-to-metals ratio similar to that of the MW was found, sug- 
gesting a higher degree of chemical evolution than in other quasar 
absorption line systems (Vladilo 2004). The extent of depletion of 
refractory elements in the absorbers with highest equivalent width 
(EW) of the Ca ii A3935 fine (1¥a3935) was shown to equal that 
seen in the warm neutral medium of the MW, in contrast to average 
DLAs. 

In this paper we present the detection of the 

[O II] AA3727,3730 emission doublet associated with Ca II 
absorbers. The detection of [O II] confirms the presence of nearby 
star formation and allows, for the first time, a direct estimate of 
the average SFR in an absorption-selected galaxy population. In 
order to provide a comparison to a more commonly studied sample 
of absorption line systems, we also analyse a sample of strong 
Mg Il-selected absorbers in which 50% are expected to be DLAs. 
Again, [O II] AA3727,3730 emission is detected. 

The outline of this paper is as follows. In Section 2 we re- 
view the sample selection, extending our original Ca ii and Mg ii 
absorber samples to lower redshifts and lower signal-to-noise ratio 
threshold. We also detail the methods used to create the composite 
spectra and measure the [O ii] emission line luminosities. In Sec- 
tion 3 we present the measured [O ii] line luminosities and discuss 
the magnitude of corrections required to convert these measure- 
ments into a total SFR for the absorber host galaxies. Wc calculate 
lower limits on the SFR per unit area of the absorbers, based simply 
on our observations and the size of the SDSS fibre aperture. Esti- 
mates of the SFR per unit area and the contribution of the absorber 
host galaxies to the volume averaged SFR density (p*) are derived 
in Section 4 which also includes a discussion of the combined un- 
certainties involved in these measurements. Finally, Section 5 dis- 
cusses the implications of the observational results for the applica- 
bility of the Schmidt law (Schmidt 1959), relating the gas surface 
density to the SFR, to the absorbers and how our results compare 
with other recent work. Unless otherwise stated, a flat cosmology 
with = 0.7, = 0.3, Ho = 100 ftkms"^ Mpc"^ and, 
where necessary, h = 0.7 is assumed throughout the paper. 



2 METHOD 

Quasar spectra were selected from the SDSS DR3 catalogue of 
Schneider et al. (2005) and supplemented with spectra in Data 
Release 4 (DR4, Adelman-McCarthy et al. 2006) main survey 
plates and spectroscopically classified as a quasar or high-redshift 
quasar^. Quasars in the sample were restricted to those with 
extinction-corrected point spread function (PSF) i-band magnitude 
(nripsf.i) brighter than 19. L The redshift range of the quasars was 
restricted to 0.41 < Zqso < 3.1; the lower limit set to ensure that 
the presence of the Mg ii AA2796, 2803 doublet could be used to 
confirm the Ca ii detections. To ensure our analysis was confined 
to relatively high quality spectra, the quasar spectra were required 
to possess ^ 3700 good pixels ("NGOOD") and spectral signal- 
to-noise ratio (SNR) ^ 6.0 in both "i" and "r" bands 

Quasar spectra showing evidence for Broad Absorption Line 
(BAL) and strong associated absorption features were removed 
from the sample. This was achieved using a preUminary Principal 

^ Component vacuum wavelength.s of 3727.09 A and 3729.88 A. 
^ 'SpecClass' equal to 3 or 4 in an SQL search of the DR4 catalogue 
■* A further ~ 150 spectra suffering from a variety of artifacts, including 
many spectra in spectroscopic plates 426 and 946 included in the DR3 re- 
lease, were removed from the sample following visual inspection. 
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Component based identification sctieme similar to tliat of WHP06, 
followed by confirmation through visual inspection. Comparison 
of our BAL and strong absorber catalogue to that of Trump et al. 
(2006) shows the two agree extremely well; it was not possible to 
simply adopt the Trump et al. (2006) DR3 catalogue because our 
sample includes quasars in the DR4 release. While it is important 
to remove the 'BAL' -quasars from our sample to obtain accurate 
dust reddening estimates, none of the results described in the paper 
are sensitive to the exact identification scheme used. In the DR3 
and DR4 samples of 27 010 and 7653 quasar spectra, 1812 and 573 
objects respectively were identified as possessing BALs or strong 
associated absorption, leaving a combined sample of 32 278 quasar 
spectra in which absorption systems were sought. 



2.1 The absorber sample 

A matched-filter search was employed to detect 
Mg IIAA2796, 2803 absorption doublets with SNR>6 over the 
redshift interval 0.4< Zghs < 1.3 in the quasar spectra®. The search 
for Mg II doublets was restricted such that the absorbers appeared 
at wavelengths > 1250 A in the quasar rest frame, i.e. not in the 
Lyman-Q forest, and had redshifts Zabs < Zqso — 0.03. The quasar 
spectra were then searched independently for Ca ii AA3935,3970 
doublet absorbers with rest frame Wxsgss > O.SA and SNR^ 
4 over the redshift interval 0.4 < Zi^hs < 1-3. The impact of 
imperfect removal of the OH skylines in the red half of the SDSS 
spectra was minimised by employing the routine of Wild & 
Hewett (2005b), with wavelength regions potentially occupied by 
known absorption lines masked. Further details of the line-finding 
procedure are given in WHP06. 

The distribution of 2:abs(Mg II)— Zabs(Ca II) shows a very 
strong excess of absorbers centred on zero- velocity difference with 
a dispersion of 70kms~^. An initial selection of 360 candidate 
Ca II absorbers was defined by selecting Ca ii absorbers within 
a velocity difference of ±200 kms^^ of a detected Mg II absorber. 
The fraction of spurious Ca ii absorbers can be estimated directly 
from the incidence of apparent Ca ii absorbers with velocity dif- 
ferences larger than ±200kms~^. The level of contamination is 
well-determined and only 5% of the sample of 360 candidate ab- 
sorbers are predicted to be due to chance coincidence. 

Absorption line properties of the Fe II A2600, 
Mg II AA2796,2803 and Mg iA2853 features for all of the 
Mg II absorber candidates were calculated via a fully parame- 
terised fit to the spectra^. The Ca II doublet was included for 
the subsample of Ca ii absorbers. A continuum was fitted to the 
regions around the Mg ii and Ca ii lines and the corresponding 
portions of the quasar spectra were normalised by dividing by 
the continuum level. Gaussian doublets were then fitted to the 
normalised spectra using a maximum-likelihood routine; the 
position and line width of the doublet, and the relative strengths 
of the individual lines, were allowed to vary freely. Rest frame 
equivalent widths (W) were calculated from the parameters of the 
Gaussian fits and errors estimated by propagation of the parameter 
errors derived during the maximum likelihood fit. 

Each candidate Ca II absorption line system was inspected vi- 
sually and 15 spurious detections were removed. Following Rao 



^ The SNR is calculated from the combined detection of both members of 
the doublet. 

^ MPFrr, http://cow.physics.wisc.edu/~craigm/idl/, Craig Markwardt IDL 
library 
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Figure 1. The wavelength region 3655 — 3795 A of a composite of the 
345 Ca II absorbers (Upper Panel) and of a composite of 3461 Mg II ab- 
sorbers (Lower Panel) in the absorber rest frame. The individual spectra 
have simply been moved to the absorber rest frame, normalised and aver- 
aged together. 



et al. (2006), Mg II systems were restricted to those with 1 < 
W^A2796/W^A2600 < 2 and M^A2796 > 0.6 and only those systems 
with Fe IIA2600 and Mg IIA2796 lines detected with a SNR>3 
(see Equation 1). This sample, referred to from now on as Mg ii- 
selected DLA candidates, was selected to maximise the fraction of 
DLA systems it contains, about 50%, thus providing a comparison 
to previous results for DLAs discussed in Section 1, and to the Ca ii 
systems which are expected to be DLAs. Our final absorber cata- 
logues contain 3461 Mg ii absorption line systems and 345 Ca ii 
absorbers in the redshift range 0.4 ^ Zabs ^ 1-3. 

The samples are more extensive than those reported in WHP06 
due to (i) the increased redshift interval searched, ^;abs > 0.4 com- 
pared to Zabs > 0.84, (ii) the larger number of quasar spectra used 
by reducing the spectrum SNR threshold from 10.0 to 6.0, and (ill) 
the lower significance threshold for the detection of the Ca II ab- 
sorbers, reduced from SNR> 5.0 to SNR> 4.0, subject to the pres- 
ence of Mg II absorption within ±200 kms~^. 

Before considering the detailed procedures employed to ob- 
tain measurements of the [O II] AA3727,3730 emission associated 
with the samples, it is worth stressing the potential of the SDSS 
spectroscopic database for such investigations. The typical la noise 
in the SDSS quasar spectra is ~ 1 x 10~^^ergscm~^ per 
pixel (Av = 69km s^^), and detection limits for (close to) un- 
resolved features in several hundred or more stacked spectra can 
easily reach 10~^*ergs cm~^ s~^. At redshifts of a~l such fluxes 
translate to unattenuated [O ii] emission corresponding to SFRs of 
significantly below 1 Moyr"^. Figure 1 shows the wavelength re- 
gion 3655 — 3795 A of composites of the Ca II and Mg Il-selected 
DLA absorbers in the absorber rest frame; the individual spectra 
have simply been moved to the absorber rest frame, normalised and 
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averaged together. The clear presence of the [O II] emission line is 
striking, including the partial resolution of the doublet. 

2.2 Building tlie composite spectra 

For the quantitative investigation of the [O ii] emission line proper- 
ties of the absorbers, the spectra were first corrected for the effects 
of Galactic extinction, then moved to the absorber rest frame, with- 
out rebinning, correcting the flux into rest-frame per-A units. Con- 
tinua in the region of [O ii] were defined using a running median 
filter (size 61 pixels) and subtracted from each spectrum. The resid- 
ual spectra were weighted by the square of the luminosity distance 
to the absorber and combined into a composite using an arithmetic 
mean ^. 

The absorber sample covers an extensive redshift interval, 
over which it is believed the average SFR in galaxies changed sig- 
nificantly (Hopkins 2004). Composites were therefore created in 
two redshift bins, Zabs = 0.4 — 0.8 and Zaha = 0.8 — 1.3. For 
the Ca II absorbers these redshift subsets were further split into 
two by Wa393b, because of the known variation of dust properties 
with strength of Ca ii absorption (Wild & Hewett 2005a). An EW- 
limit of Wa3935 = 0.68 A was chosen for this division to match 
that used in WHP06 but the results reported below are not sensitive 
to the exact value of EW used. Unlike any of the Ca ii absorbers, 
~ 0.2% of the Mg II absorber spectra show detectable [O ii] emis- 
sion in individual spectra. However, the contribution to the mean 
[O II] luminosity of the absorber sample by these objects is small 
and the results of the paper are not dependent on the inclusion of 
the small percentage of such systems. The composite spectra in the 
region of the [O ii] line are shown in Figs. 2 and 3. 



tween the samples split by redshift, particularly for the Ca II ab- 
sorbers, which may be caused, in part at least, by aperture bias 
— less galaxy light falling within the fibre radius at low redshifts. 
Secondly, there appears to be a trend of increasing SFR with in- 
creasing Wa3935 in the low-redshift Ca ii subsample, but no trend 
is detected at high redshift, although the errors are large. Again, 
if stronger Ca II absorbers were more likely to be found at smaller 
impact parameters from the galaxy centre, the difference could be a 
manifestation of aperture bias. Due to the large errors and evidently 
weak, or non-existent, trends with equivalent width, the remainder 
of this paper will focus on the samples split by redshift only. 



3 STAR FORMATION RATES 

We convert the [O II] line luminosities into instantaneous SFRs, 
listed in column 6 of Table 1, using the relation of Kewley et al. 
(2004): 

SFR([On](M0yr-^) = 6.58 ± 1.65 x 10"*^i([On]) (2) 

where [O II] luminosity is measured in ergs s^^. The errors on our 
measured line fluxes and the uncertainty associated with the SFR 
conversion have been propagated in the usual way. In the follow- 
ing subsections we discuss three correction factors which must be 
applied to these values in order to estimate the true average SFR of 
absorber host galaxies. The first two are caused by dust within the 
ISM of the absorber host galaxy, the third by the finite aperture of 
the SDSS fibres. 



2.3 [O II] AA3727,3730 line fitting 

In each composite spectrum the [O ii] AA3727,3730 doublet was 
fitted with a double Gaussian with line ratio fixed at 4:5®, using a 
non-Unear least squares fitting routine, and line limiinosities mea- 
sured from the fitted Gaussian parameters. Whether or not a fea- 
ture was detected was determined via its SNR (Hewett et al. 1985; 
Bolton et al. 2004, for a full derivation): 



SNR = 



(1) 



where is the composite spectrum and a is the flux error, w is a 
Gaussian line profile of position and width given by the fitted line 
parameters. A positive detection was defined to have a SNR> 3.0. 
Where a line was undetected, upper limits were placed for a single 
Gaussian line of a twice the pixel resolution and detection limit of 
SNR= 3.0. 

Table 1 presents the measured [O ii] line luminosities and er- 
rors for each composite. From these values we make a few im- 
mediate observations. Firstly, there is a significant difference be- 



In this situation it is not possible to use a weighted mean due to the cor- 
relation between absorber redshift and SNR of the spectra: low-redshift 
absorbers can be detected in lower-redshift quasars which have generally 
higher SNR spectra. 

The 4:5 line ratio is that predicted for T=10 0()0K and electron density 
A^e = 200 cm^^ (Osterbrock 1989) and it provides a good fit to the com- 
posite [O II] AA3727,3730 doublet. The line strength measures are not sen- 
sitive to reasonable variations in the adopted line ratio. 



3.1 Dust corrections to Ca ii absorbers 

While the dust content of the Mg ii-selected DLA candidates is 
negligible (WHP06), it is necessary to correct our measured SFRs 
for the effect of dust in the Ca ii host galaxies. This correction is 
twofold. Firstly, the [O II] photons are absorbed by dust grains, thus 
reducing the line flux emitted from the galaxy. Secondly, the dust 
obscures the background quasars, thus reducing the total redshift 
path available to find Ca II absorption line systems in the survey 
over that naively expected. To estimate the magnitude of both ef- 
fects we require the extinction as a function of wavelength caused 
by the dust. This can be achieved from the measured reddening in 
each subset, by assuming the form of the extinction curve and the 
total-to-selective extinction ratio (B.v = Av /E(B—V)). In what 
follows we assume the dust content as measured from the absorp- 
tion line-of-sight is applicable to that obscuring the nebular emis- 
sion line region. In Section 4.4.2 we discuss the potential problems 
with this assumption given the likely complicated geometry of the 
dust distribution within galaxies. 

For a full description of our method to measure dust red- 
dening in the absorbers see WHP06; a brief description is given 
here. High SNR quasar template spectra are created by combining 
the spectra of all quasars in our input sample which do not con- 
tain known Mg II systems, a total of 22 727 quasar spectra with 
0.35 < Zqso < 3.15. The template spectra are created in redshift 
bins of A^qso = 0.1. Each quasar spectrum containing an inter- 
vening absorption line system is then divided by a suitably nor- 
malised high SNR quasar template spectrum. Each residual spec- 
trum is then moved to the rest frame of the absorber and all residual 
spectra in the subsample are combined using an arithmetic mean. 
The resulting residual spectrum is fitted with the Large Magellanic 
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Figure 2. SDSS composite spectra of Ca II absorption line systems (top) and Mg Il-selected DLA candidates (bottom) in the region of [O II] AA3727,3V30 
in redshift bins of 0.4 z^bs < 0-8 (left) and 0.8 ^ Zabs < 1-3 (right). Flux from each contributing spectrum has been converted into luminosity before 
combining. Overplotted (in red lines) are Gaussian fits to the emission lines. Measured fine luminosities from these line fits are given in each plot. 
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Figure 3. Same as Fig. 2 for the four subsamples of Ca II absorbers split by redshift and W;^3935. No line is formally detected in the top left plot (SNR< 3). 
Note the change in y-axis scale between the upper and lower plots. 
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Table 1. Sample sizes, reddening, measured [O II] line luminosities, inferred SFRs before and after correction for dust, and dust obscuration bias 
effecting the observed population, for each Ca II composite spectrum and the two Mg ll-selected DLA candidate composite spectra. 



Sample 


Nabs 


E(B-V) (mag) 


L{0 II] (10*0 erg s-i) 


SNR 


SFR(M0yr-l) 


SFR" (]V[0yr-l) 


Dust bias'' (E(B-V)ma.x) 


Low z 


215 


0.036 ± 0.005 


1.4 ±0.2 


5.9 


0.09 ± 0.01 


0.11 ±0.02 


0.24 (0.2) 


High 2 


130 


0.049 ± 0.004 


5.8 ±0.6 


5.8 


0.38 ± 0.04 


0.48 ± 0.05 


0.45 (0.25) 


Low z. Low EW 


153 


0.029 ± 0.006 


< 0.9 


3.0 


< 0.06 


< 0.07 


0.20 (0.2) 


Low 2, High EW 


62 


0.051 ± 0.009 


2.7 ±0.4 


6.6 


0.18 ±0.03 


0.22 ± 0.03 


0.33 (0.2) 


High z. Low EW 


74 


0.036 ± 0.006 


5.4 ± 1.3 


4.3 


0.35 ± 0.09 


0.41 ±0.10 


0.34 (0.25) 


High z. High EW 


56 


0.066 ± 0.007 


6.6 ± 1.2 


5.0 


0.43 ± 0.08 


0.58 ±0.10 


0.56 (0.25) 


Low z, Mg II 


1291 




1.7 ±0.06 


15.2 


0.11 ± 0.004 






Higli Mg II 


2170 




2.1 = 0.:-! 


7.0 


0.11 ±0.02 







" Corrected for dust attenuation at [O II], based on E(B—V) derived from reddening of the background quasar. 

^ Fraction of the total population of Ca II absorbers lost due to dust obscuration of the backgrovmd quasar sample. This calculation is only performed 
for absorbers with E(B—V) less than the value given in brackets. 



Cloud (LIVIC) reddening curve of Pel (1992) and the resulting 
E(B—V) = Ab — Av values are given in the third column of Ta- 
ble 1. 

The formal errors on tMs fit are small and the dominant source 
of error arises from the intrinsic variation in quasar spectral en- 
ergy distributions (SEDs). This error is estimated via 10 000 Monte 
Carlo simulations drawing a random subset of quasar spectra equal 
in size to the subset of absorbers. A random absorber redshift is 
assigned to each quasar spectrum and the above method repeated 
to estimate the reddening of this simulated sample of absorbers. 
Repeating 10000 times, the 68th percentiles of the final reddening 
distribution are taken to provide the errors quoted in the third col- 
unm of Table 1. Again, further details on this error calculation are 
provided in WHP06. 



3.1.1 [Oii\ line attenuation 

The attenuation of light due to dust at the wavelength of [O II] 
is calculated from the measured E(B—V) values in the absorption 
line systems, and the assumed extinction curve. The resulting dust 
attenuation corrected [O ii] SFRs are given in colunm 7 of Table 
1. Again, the errors have been propagated in the usual way. In the 
calculation we have assumed that the photons from the star forming 
regions are passing through the same quantity and type of ISIVI as 
the light from the background quasar, an assumption that will not 
be correct in detail. However, without further information on the 
nature of the relation between the absorbers and host galaxies it is 
difficult to improve on this assumption. The effect on our results of 
possible enhanced extinction in the nebular emission line regions is 
discussed further in Section 4.4.2. 



3.1.2 Dust obscuration of the background quasars 

The second effect of dust in the absorber is to obscure the back- 
groimd quasars, meaning that fewer Ca ii absorption Une systems 
are found in a magnitude Umited quasar sample than if the ab- 
sorbers did not contain dust. The method of correction for this "dust 



® WHP06 found the LMC reddening curve provided the best fit to the com- 
bined quasar residual spectra, although the derived values of E(B—V) are 
not sensitive to the exact choice of extinction curve. 



obscuration bias" is discussed in detail in WHP06, but given the im- 
portance of the correction we include a revised description of the 
procedure in Appendix A. 

The fraction, (6), of the total Ca ii absorber population lost 
due to dust obscuration bias is given in the final column of Table 1 
for each Ca II sample. To obtain the true number of Ca II absorbers 
which would have been observed had they contained no dust, the 
observed number is divided by (1 — 6). The redshift, as well as the 
dust content, of the absorbers has a large effect on the extent of the 
bias due to the shape of the dust extinction curves — the strongest 
attenuation is at small wavelengths; higher redshift absorbers there- 
fore cause greater extinction to the background quasar. The smallest 
correction factor applies to the least dusty absorbers at the lowest 
redshift. Conversely for the absorbers with Wx^nr, > 0.68 and 
0.8 S5 z^bs < 1.3 we are in fact only observing 44% of the to- 
tal number of Ca II absorbers with E(B—V) < 0.25. As argued in 
WHP06, these correction factors are very well determined. How- 
ever, as discussed in Appendix A, there is an important caveat; we 
are unable to correct reliably for absorbers with an E(B—V)> 0.2 
and > 0.25 for the low- and high-redshift samples respectively be- 
cause our observations have little or no sensitivity to absorbers with 
such values of E(B—V). 



3.2 Non-DLA contamination of IVIg ii-selected DLA sample 

Bias due both to [O ii] line attenuation and the obscuration of back- 
ground quasars is minimal for the Mg ii-selected DLA candidate 
sample given the very low dust content (see WHP06) and we have 
not applied any corrections to the measured line fluxes. One un- 
certainty in the implied SFR of DLAs from our analysis lies in 
the contribution made by the 50% non-DLA contaminants to this 
sample. An estimate of the SFR of these contaminants can be cal- 
culated by measuring the SFR found in IVIg II absorbers which fall 
out of the predicted DLA Fe ii/JVIg II EW ratio range, but still with 
Wa2796 > O.GA. In these absorbers we measure SFRs of 31(57)% 
those found in the low- (high-) redshift Mg ii-selected DLA can- 
didate sample, producing corrections to the SFRs quoted in the fi- 
nal two rows of Table 1 of a factor of 1.69(1.43). We note that this 
method of estimating the correction factors is only approximate and 
the associated uncertainties are discussed further in Section 4.4. 
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Figure 4. Proper transverse distance covered by a 3 arcsec diameter fibre 
as a function of redshift, for the ACDM cosmology used in tliis paper with 
Ho = lOOkms-iMpc-i. 



3.3 Aperture Corrections 

The final correction factor which must be applied to our measured 
SFRs in both the Ca II and Mg II samples is for the finite aperture 
of the SDSS fibre, which is centred on the background quasar and 
not on the host galaxy of the absorber. The fibre may be covering a 
fraction of empty sky, which has implications for our measures of 
SFR per unit area of the absorbers; or it may be missing a fraction 
of the galaxy light, with implications for our estimate of the SFR 
per absorber; or both. 

For ease of reference, Fig. 4 depicts the proper transverse dis- 
tance covered by a 3 arcsec diameter fibre as a function of red- 
shift. Atmospheric seeing will cause a blurring of the discontinuity 
at the edge of the fibre, but in practice star formation occurring 
within ~ 17/2 = 8.5 /i~^kpc from the line-of-sight to the quasars 
will contribute to the observed SFR for our highest redshift ab- 
sorbers. At the lowest redshifts in our sample this is reduced to 
~ 14/2 = 7 /i~^kpc. The fibres cover a relatively large area, com- 
parable to the size of the star forming extent of a galaxy at these 
redshifts (Ferguson et al. 2004), but a smaller area than the likely 
extent of, for example, Mg II or H I halos/disks surrounding the 
galaxies. To calculate many quantities of interest, such as SFR den- 
sity per unit volume as a function of redshift, p* (z), information on 
the size of the absorbers is required and this aspect of the discussion 
is deferred to Section 4. 



3.4 Summary of observational results 

Table 1 summarises the principal observational result of this paper 
relating to the integrated SFR for the absorbers, i.e. the determina- 
tion of the observed average SFR for the Ca ii absorbers and the 
Mg ii-selected DLA candidates. In the remainder of this section 
we look at the SFR per unit area of the absorber host galaxies, for 
which a hard lower limit can be obtained by dividing the observed 
average SFR per absorber by the area of the spectroscopic aperture 
at the absorber redshift. The lower limit arises due to the uncer- 
tainty as to the fraction of the fibre which falls on blank sky. 



3.4.1 Mg u-selected DLAs 

By dividing the SFRs given in the final two rows of Table 1 by the 
proper area of the fibre at the central redshift of the samples, and 
correcting for the 50% of non-DLAs in the sample as discussed 
above, we obtain limits on the SFR per unit area of DLAs. The lim- 
its are > 12 ± 0.4 X 10"* and > 9 ± 1 X lO""* ft^ Mgyr-^ kpc"^ 
for the low-z ({z) ~ 0.6) and high-z ({z) ~ 1.0) samples. 

To put these limits into context, we compare with the Kenni- 
cutt parameterisation of the Schmidt law (Schmidt 1959; Kennicutt 
1998a; 1998b) which relates the surface densities of neutral gas and 
SFR in galaxies: 



^ = 2.5 X 10"* X ( 



iV(HI) 



SFR 



1.25 X 1020 cm-2 



(3) 



above a critical H I column density, Nc = 5 x 10^" cm"^^ where 
XisFR u"its of Moyr~^ kpc~^. A distribution function for 
the H I column densities of the absorbers is thus required to pre- 
dict their SFR per unit area. In fact, the combination of the critical 
threshold value Nc for the onset of star formation, combined with 
the non-linear dependence of the SFR on column density, means 
that simply evaluating Equation 3 for (A'^(H I)) produces an answer 
accurate to better than 10%. 

Taking the distribution of H I column densities for DLAs at 
1.7 < z < 2.2 of Prochaska et al. (2005) results in a mean column 
density of 6.6 x 10'^'^ cm"^. The mean is insensitive both to the ex- 
act form of the acceptable functional fits given in Prochaska et al. 
(2005) and the upper cutoff in H I column density between 10^^ and 
lO'^^ cm~2. The latter behaviour is due to the steep decline in the 
number of systems at very high column densities. Unfortunately, 
at « < 2, the redshift range applicable to our systems, significant 
imcertainty in the predicted SFR per unit area for the absorber pop- 
ulations comes about because of the poorly determined distribution 
of iV(H I). Rao et al. (2006) find that the mean iV(H I) is somewhat 
higher than assumed above, (iV(H I)) = 1.27 ± 0.36 x 10^^ cm'^ 
for DLAs with O.IK « < 0.9, and 1.07 ± 0.23 x 10^^ cm'^ for 
DLAs with 0.9 < z < 1.65. 

The Rao et al. (2006) mean column density (N(ll I)) cor- 
responds to a SFR per unit area of 13.1 ± 5.2 x 10""^ and 
10.3 ± 3.1 X 10"^ /i^ Moyr"^ kpc"^ according to Equation 3, 
for their low and high redshift bins respectively. These values are 
factors of 1 1 ± 4 and 1 1 ± 3 higher than the lower limits we deduced 
for our low- and high-z samples (the errors here are dominated 
by the uncertainty in {NQi I))). For comparison, the high redshift 
(A'^(H I)) value of Prochaska et al. (2005) results in a predicted 
SFR per unit area for DLAs of 5.2 x 10"^ h'^ Mgyr"^ kpc"^. 
This is still 4(6) times greater than the derived lower limit in the 
low- (high-) z Mg Il-selected DLAs. 

Does the mean H I column density of DLAs apply to our 
Mg ii-selected sample? These absorbers were selected to maximise 
the DLA fraction they contain and allowance was made to correct 
the observed SFR for contamination by non-DLAs (Section 3.2). 
The results of Rao et al. (2006) are based on follow-up A'^(H I) ob- 
servations of a sample of Mg ii-selected systems very similar to 
that used here and it would be very surprising if the properties of 
our Mg II absorber sample were significantly different from those 
studied by Rao et al. (2006). 



3.4.2 Ca II absorbers 

In the case of the Ca ii absorbers, dividing the SFRs by the proper 
area of the fibre at the central redshift of the samples results in lim- 
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its for the SFR per unit area of > 7(21) x 10"" Mc5yr"^ kpc"^ 
for the low- (high-) z samples respectively, after correcting the 
SFRs for attenuation of the [O ii] line by dust. 

It is difficult to compare these results to those predicted by the 
Schmidt law, due to the very few systems currently with both Ca ii 
and A'^(H I) measures, particularly for the high IFAaoar) systems. 
While the Ca ii absorbers are expected to have N(]l I) values above 
the nominal Umit for DLAs, the comparison to predicted SFRs per 
unit area of DLAs should be treated with caution. 

In the following section, we discuss the magnitude of correc- 
tions applicable to these values caused by portions of the fibre aper- 
ture covering blank sky, and convert the measured lower limits into 
absolute SFR per unit area values for the absorbers. Here we sim- 
ply note that the corrections are expected to be small, particularly 
in the case of the Mg ii-selected DLAs. 



4 ABSORBER STAR FORMATION PROPERTIES 

In this section we use results of imaging surveys of DLAs to esti- 
mate their average extent and thus assess the magnitude of correc- 
tions applicable to our observed SFRs within the fibre aperture. The 
discussion in this section is predicated on the assumption that the 
SFR per unit area is determined entirely by the Schmidt-law via the 
H I column density, as specified in Equation 3. If DLAs at 2; < 1 
are caused by large gas disks surrounding a central galaxy, the ab- 
solute SFR per unit area observed through the SDSS fibre aperture 
can be predicted from the mean N{Yi 1) of the population and, fur- 
thermore, we can estimate the contribution of the host galaxies to 
the star formation rate density (p*) of the universe. Discussion of 
scenarios in which the relative location of star formation and ab- 
sorption within a galaxy/halo system, and/or failure of the Schmidt 
law, impact greatly on our results follows in Section 5. 

4.1 Physical size of absorbers 

In the case of DLAs, imaging of host galaxies at « < 1 suggests 
they may occur in halos or disks that extend out as far as ~ 25/i~^ 
kpc from the host galaxy (e.g. Chen & Lanzetta 2003; Rao et al. 
2003; sec Tabic 2 of Zwaan et al. 2005 for a recent compilation). 
The mean observed impact parameter between DLAs and absorber 
host galaxies is ~ 9h~^ kpc (Table 2 of Zwaan et al. 2005). Assum- 
ing a simple circular geometry for the absorbers this corresponds 
to a typical absorber radius of ~ 13/i~^ kpc, assuming covering 
factors of unity. Comparatively, from direct H I observations of lo- 
cal galaxies, Zwaan et al. (2005) predict a very similar mean impact 
parameter of ~ 7h~^ kpc. 

We therefore adopt 13/i~^kpc as the predicted extent from 
the centre of the host galaxy of absorbing clouds which give rise to 
DLAs. Adopting the measured value of dP/dz = 0.079 ± 0.019 
(Rao et al. 2006) for DLA absorbers at low redshifts, the estimated 
size leads to absorber space densities of 0.027(0.022) Mpc~^ 
for our low- and (high-) redshift samples (Equation 5). 

We turn now to the Ca 11 absorbers. Without any direct obser- 
vations of host galaxies, or information on how far Ca 11 extends 
into the gaseous halos surroimding galaxies, we must make some 
assumption about their size. In WHP06 strong arguments are given 
that Ca II absorbers have N(H i) values above the nominal limit 

The mean observed impact parameter is about two-thirds of the maxi- 
mum radius, assuming the absorbers present simple circular cross-sections. 



for DLAs and. as such, make up a subsample of DLAs. Their high 
dust and metal contents suggest a scenario in which Ca 11 absorbers 
exist closer to the central galaxy, possibly with a significant frac- 
tion lying within the star forming region of the galaxy. The values 
of dP/dz {Wx379s > 0.5 A) for the Ca 11 absorbers, correcting 
for those missed from the sample duo to dust obscuration bias, arc 
0.019 and 0.025 for our low- and high-redshift samples (WHP06 
and see Section 4.3.1). Taking the same space densities as for the 
DLA absorbers and simply scaling the cross-sections according to 
the difference in the observed dP/dz for the two absorber popu- 
lations gives radii for Ca II absorbers of 0.49(0.56) times those 
of DLAs, or 6.4(7.3)/i"^kpc for the low- (high-) z samples". We 
note that the difference in size for the two redshift intervals is not 
significant given the combined errors on number density for both 
the DLAs and Ca 11 absorbers. 



4.2 SFR per unit area 

Given the physical size of the absorbers we can calculate their ac- 
tual SFR per unit area, within our assumed geometrical model, for 
which lower limits were calculated in Section 3.4. 

Beginning with the Mg ii-selected DLA candidates, plac- 
ing 6 — 9/i~^kpc radius circles (i.e. the SDSS spectroscopic fi- 
bres, see Section 3.3) at random across a second circle of ra- 
dius 13/t"^kpc (i.e. the gaseous halos or disks), such that the im- 
pact parameter between the centre of the circles never exceeds 
the radius of the second circle, and simmiing the fraction of area 
falling outside the absorber boundary (i.e. on blank sky) leads to 
aperture corrections of 1.16 and 1.21 for low- and high-z sam- 
ples respectively. Thus, the SFR per unit area of Mg ll-selected 
DLA candidates derived in Section 3.4 can be increased only to 
14±0.5(11 ± 1.6) X 10"* Mgyr"^ kpc"^ and alarge discrep- 
ancy between the conventional Schmidt-relation and our observa- 
tions remains. 

Adopting the estimates for the size of the Ca II absorbers 
leads to an aperture correction of 1.59 (1.69) for the low- (high) 
z samples. These corrections lead to a SFR per unit area of 
11(36) X 10"'' Moyr"' kpc"^, i.e. very similar to the Mg II- 
selected DLA candidates at low redshift, but four times higher in 
the high redshift sample, //'the mean H I column of Ca II absorbers 
is equal to that of DLAs with 1.7 < z < 2.2, the high redshift 
Ca II absorber host galaxies have SFRs per unit area closer to that 
predicted by the Schmidt law, while at low redshifts an order of 
magnitude discrepancy remains. 



4.3 The global SFR density 

The global SFR density of the universe has now been measured by 
many surveys out to high redshift. The general consensus is that, 
as we look back in time, the SFR density increases by about one 
order of magnitude from z = to 1, and remains approximately 
constant at higher redshifts up to 2: ~ 5 (Bouwens & lUingworth 
2006). Hopkins (2004) has provided a power-law fit to the data at 
z<l: 

log p* = 3.29 log(l + z)-1.8 (4) 



Despite the different method used to estimate the halo sizes of Ca II ab- 
sorbers, the resulting size remains very similar to that suggested in WHP06. 
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where p* is in units of MQyr^^ Mpc^'' and Ho = 
70 kms~^ Mpc^^. In our low- (high-) z bins the average SFR den- 
sity of the universe is thus log p* ~ —1.13 (—0.81) or, with p* in 
units of M©yr"^ Mpc"^, log p* ~ -0.66(-0.35). From now 
on all p* values will be quoted in the latter units. 

Adopting either a physical size or a space density for the Ca II 
and Mg ii-selected DLA absorbers, and assuming the star forma- 
tion is uniformly distributed (on scales similar to the fibre size) 
within the high column density gas, as required by the Schmidt 
law, we can estimate the contribution to the global SFR density of 
the universe from the host galaxies of each class of absorbers. 

4.3.1 Cross sectional area and comoving space density of the 
absorbers 

The observed redshift number density (dP/dz) of absorbers pro- 
vides a joint constraint on their size and space density. For ab- 
sorbers with cross-sectional area, S, and comoving number den- 
sity, n{z), their number per unit redshift is given by (Hogg 1999; 
Peebles 1993): 

-;— = S n(z) — , (5) 

^ ' Ho ^Um{1 + zY + Q.K 

The dP/dz of the whole Ca ii absorber population, i.e. cor- 
recting for dust obscuration bias (Section 3.1.2), is obtained by 
multiplying the observed dP/dz by the fraction lost due to dust 
bias (final column of Table 1). For the low- and high-redshift sam- 
ples we find dP/dz (M^x3795 > 0.5 A) = 0.019 and 0.025 respec- 
tively^^ (see WHP06 for a detailed method). In Fig. 5 we plot the 
combined constraint on comoving number density and cross sec- 
tional area of the dust bias corrected Ca II absorber samples, given 
their respective redshift number densities. The constraints for two 
different redshifts within each redshift interval are plotted, although 
the change within a bin is not large. 

To provide reference points, the z^Q space density of galaxies 
with S-band magnitudes, M > M* + 1.0 and M > M* + 2.0 
are indicated by horizontal lines. M* = —19.79 + 51ogio(ft) 
and 4>* = 0.00159/1^ Mpc~^, the characteristic magnitude (on 
the Vega-system) and normalisation of the luminosity function, are 
taken from the rest-frame B-band Schechter function parameterisa- 
tion of the 2dF Galaxy Redshift Survey (2dFGRS, Madgwick et al. 
2002). 

4.3.2 Contribution of absorber host galaxies to p* 

Starting with the Ca ii absorbers, within the errors, the space den- 
sity estimates of 0.027(0.022) /tnipc"'' (Section 4.1), for the 
low- (high-) z samples, are consistent with zero evolution and are 
similar to the space density of galaxies brighter than Ml^p -|- 2. 

Given the small impact parameters expected for the Ca ii ab- 
sorbers within our assumed model for the relative configuration 
of star formation and absorbing gas, the central host galaxy is 
expected to fall within the fibre aperture essentially all the time 
(Fig. 4). The fibre aperture correction factor (that quantifies the 
host galaxy emission line flux lying outside the fibre diameter) is 
1.32(1.25) for the low- (high-) z samples. The volume averaged 

^'^ Strictly speaking these are lower limits to the true number densities, as 
we are unable to correct for dust obscuration bias above the iJfi?— Vj^^^, 
values given in Table 1. 



SFR density of the absorber host galaxies, obtained from the prod- 
uct of the SFR, the space density and the aperture corrections, is 
then log p* ~ — 2.43(— 1.89). Comparing to Equation 4 these val- 
ues represent only 2(3)% of the total SFR density in the universe at 
the same redshifts. 

Turning now to the Mg Il-selected DLA absorbers, dP/dz 
is a factor three to four larger and, unless the DLAs possess a 
much higher space density than that of the Ca ii absorbers, allow- 
ing them to be considerably smaller, their physical extents are such 
that a significant aperture correction is necessary in order to al- 
low for the absorber gas that falls outside the 3 arcsec SDSS fibre 
aperture. Adopting the absorber size estimate of 13h~^ kpc, and 
hence space density, from Section 4.1 leads to aperture corrections 
of 3.9(2.9) for the low- (high-) z samples respectively. Taking the 
measured SFR for the Mg ii-selected DLA candidate sample then 
leads to values of logp* — — 1.72(— 1.91). These values corre- 
spond to ~9(3)% of the SFR density in the universe at <z>~0.6 
and <a>~1.0 respectively. 

It is not possible to explain the small contributions to the 
global SFR density by simply appealing to an underestimate of the 
size of the absorbers in Section 4.1. Once absorber size is compa- 
rable to or larger than the fibre-size (as in all examples here), and 
assuming the Schmidt law holds, changes in the assumed absorber 
size make no significant difference to the value of p* . This follows 
from p* being proportional to the product of the comoving number 
density and aperture correction: 

p* oc (SFR) X n{z) x A (6) 
oc (SFR)x^x^ 

where (SFR) is the SFR per absorber, A is the aperture correction 
factor, n(z) the comoving number density, S the area of the ab- 
sorber and dP/dz is fixed by observation. In the limit that the size 
of the absorber (and star formation region) is greater than the fibre 
aperture, as S increases the aperture correction increases propor- 
tionately and p* remains constant. 

Our results for the p* of Mg ii-selected DLA candidates are 
directly comparable to those of Hopkins et al. (2005) who con- 
cluded that DLAs are responsible for as much as 80% of the SFR 
density of the universe at similar redshifts to those considered here. 
The dramatic discrepancy between our study and theirs arises di- 
rectly from the assumption by Hopkins et al. that the local Schmidt 
relation holds for DLA absorbers. Assuming the mean A'^(H I) for 
our Mg ii-selected DLAs from Rao et al. (2006), as appropriate for 
a direct comparison between the results, we observe a S¥9. per unit 
area of <10% of that expected were the Schmidt law to apply. 



4.4 Summary of uncertainties 

For each of our absorber samples potential sources of uncertainty 
in our analysis have been made clear throughout the paper The fi- 
nal SFR measurements presented in Table 1 include associated line 
measurement errors, the error on the conversion from [O ii] lumi- 
nosity to S¥R. and, in the case of the Ca ii absorbers, the error on 
the dust attenuation correction estimated from the reddening of the 
background quasar continuum. In this subsection we combine these 
errors with additional uncertainties in our calculations to produce 
an estimate of the reliability of our key results. 
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Figure 5. Combined constraints on cross sectional area and comoving number density, given the measured dP/dz of the Ca II absorption line systems 
corrected for those missed from the sample due to dust obscuration bias. Left: low redshift sample: yighr. high redshift sample. In each plot, the constraints are 
shown for two different redshifts within the redshift bin. The horizontal lines are comoving number densities of galaxies, at z ~ 0, with M > M* + 1 and 
M > M* + 2, where M" is the 2dFGRS characteristic magnitude (see text). 



4. 4. 1 Star formation rates from [Oil] emission 

Kewley et al. (2004) present a comprehensive analysis of the de- 
pendence of the [O ii]-derived SFR in systems with an extended 
range of physical properties, including objects with SFRs signifi- 
cantly lower than those we measure for the absorber sample. Given 
the known low metallicity of DLA gas, the adopted Kewley et al. 
(2004) relation applies for systems with a mean metallicity that 
may be higher than that of the absorber populations. Applying 
a correction for metallicity would decrease the inferred SFR for 
a given [O ii]-luminosity, increasing the disagreement between 
our observations and predictions based on the application of the 
Schmidt law. As very few absorption systems have measured emis- 
sion line abundances, and the link between absorption and emission 
line metallicities therefore remains unclear, it is not possible to as- 
sess this effect further at this stage. 

4.4.2 Dust attenuation of the [O ll] emission 

In calculating the SFR of the Mg ii-selected DLAs we have not 

made any correction for dust attenuation of the [O II] line, be- 
cause of the very low reddening detected in Mg ii absorbers by 
WHP06. Specifically, for an equivalently defined sample we found 
an E(B—V) of 0.008 ± 0.0016 which results in a flux correction 
of <4% at 3729A, insignificant compared to the line measurement 
errors. 

Unlike the Mg 11 absorbers, the Ca 11 absorbers do contain 
significant quantities of dust. However, as argued in WHP06, the 
quantity of dust, associated errors and resulting obscuration bias are 
well determined. The E(B—V) values obtained are not dissimilar to 
those expected from sightlines passing through the disks of spiral 
galaxies. 

However, it is possible, and perhaps likely, that nebular emis- 
sion line regions suffer significantly greater dust extinction than 
that experienced by the majority of the interstellar medium. Bell & 
Kennicutt (2001) find a median E(B~V) of ~ 0.4 for nebular emis- 
sion regions of nearby star forming galaxies, although the mean 
SFR of their sample is considerably larger than for the systems 
studied in this paper. 

With the currently available data it is only possible to obtain 
a rough estimate of the true attenuation in the nebular emission 



regions we observe. By combining all 1559 Mg II absorbers with 
Zahs < 0.85, in which H0 is potentially visible, we find an [O II] 
to H/3 ratio of 1.9. Given the theoretical balmer decrement (Oster- 
brock 1989), relation between Ha line strength and SFR (Keimi- 
cutt 1998), and between [0 11] line strength and SFR (Kewley et al. 
2004), we predict a ratio of 3.44. This suggests dust attenuation 
may cause us to be underestimating the SFRs by up to a factor of 
two, however, we emphasise that an accurate analysis must await 
further observations. 



4.4.3 Mg ll-selected DLAs 

Firstly, we look at the result with perhaps the farthest reaching im- 
plications, that the SFR per unit area in DLA gas is observed to be 
a magnitude lower than that expected from the Schmidt law. We 
consider in turn the effects of the adopted mean A'^(H I), DLA frac- 
tion, absorber size, contamination by non-DLA absorbers to the 
measured SFR and dust. 

As stated in Section 3.4 it would be very surprising were the 
mean A'^(H I) and contamination fraction of our Mg Il-selected 
DLA sample to turn out to be very different from those found by 
Rao et al. (2006) given the essentially identical criteria used to de- 
fine both absorber samples. Were we to overestimate the size of the 
absorbers we would underestimate the SFR per unit area. The sce- 
nario in which the absorber size is reduced to bring our result into 
line with the Schmidt law is discussed in Section 5 and found to be 
highly implausible. 

Turning to the contamination of the sample by non-DLA ab- 
sorbers, we concluded in Section 3.2 that non-DLAs could have 
SFRs of ~ 31 - 57% that of the DLAs. If, however, the non- 
DLAs contributed zero SFR to the measured lines our SFR per 
absorber would increase by only a further 20% and 40% for the 
low- and high-redshift samples respectively. Combining the final 
SFR per unit area of Section 4.2 with the possibility that the non- 
DLAs in our sample contribute zero flux to the observed lines 
allows a SFR per unit area of at most 16.6 ± 0.6(15 ± 2) x 
10~* M0yr~^ kpc~^, i.e. still a factor of eight (seven) below 
the Schmidt law for the low- (high-) 2 samples respectively. 

Finally, as discussed in the previous subsection, attenuation 
of the [O II] line by dust may be our biggest cause of uncertainty. 
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due to the unknown distribution of dust around star forming regions 
compared to in the interstellar medium. Allowing for an underes- 
timate of the SFR by a factor of two, and that the non-DLAs con- 
tribute zero flux, would reduce the discrepancy between our results 
and the Schmidt law to a factor of four (three). 

Turning now to the SFR density measurements, assuming the 
Schmidt law holds and allowing for the non-DLA contaminants to 
contribute zero flux to the measured lines, we find maximum pos- 
sible contributions of DLA host galaxies to the overall SFR density 
of the Universe of 10 ± 3% (4 ± 1%). Within the adopted model, 
an underestimate of the absorber size does not change our results, 
due to the compensation by the aperture correction (Equation 6). 
AUowing for the possibility that we have underestimated dust at- 
tenuation by at most a factor of two brings this result to 20% (8%). 

4.4.4 Ca II absorbers 

As with the Mg Il-absorbers, an underestimate of the absorber size 
can not change our results, provided the Schmidt law applies. Even 
given the imcertainty in dust extinction, it is difficult to see how the 
inferred contribution of Ca ii absorber host galaxy light, that falls 
within the ~ Ihh^^ kpc diameter covered by the SDSS-fibres, can 
contribute more than a few percent to the global SFR density. 



5 DISCUSSION 

Given our lack of knowledge concerning the size distribution and 
geometry of Ca II and DLA absorbers we have, in the previous sec- 
tion, deliberately confined ourselves to providing illustrative cal- 
culations of the SFR per unit area and the volume averaged SFR 
density based on conventional assumptions concerning the size of 
the absorbers and the relationship between A'"(H I) and SFR. In 
particular, no account has been taken of more realistic geometries 
that may be involved in the relation between the measured NQi i)- 
column along a single line-of-sight through an absorber and the 
spatial distribution of star formation within the system giving rise 
to the absorber. Further quantitative progress must await the results 
of imaging investigations capable of constraining the extent of the 
star forming regions in Ca ii and DLA absorbers at z<l. How- 
ever, the detection of weak [O ii] AA3727,3730 nebular emission 
associated with the Ca II and Mg Il-selected DLA candidate ab- 
sorbers places strong constraints on the combination of absorber 
size, space density and star formation properties of absorption se- 
lected systems. 

The lower value of dP/dz and the significantly higher SFR 
seen in the high redshift Ca II absorbers compared to the Mg II- 
selected DLA candidate absorbers fits in well with the view ad- 
vanced in WHP06 that the Ca ii absorbers represent chemically 
evolved systems with detectable quantities of dust that may be spa- 
tially closely related to intermediate brightness late-type galaxies. 
That the combination of absorber size and SFR produces a SFR per 
tmit area closer to that expected from the Schmidt law, assuming 
Ca II absorbers have mean iV(H i) similar to that found in DLAs, 
adds further weight to such an interpretation. However, the con- 
sequence of this picture is that Ca ii absorbers then account for 
only a small fraction, just ~ 3%, of the total SFR density in the 
tmiverse at 2 ~ 0.5 — 1. Appealing to a missing fraction of ab- 
sorbers with associated extinction exceeding our sensitivity limit 
of E(B—V) ~ 0.25 mag can increase the fractional contribution to 
the p* but not to the extent of explaining the factor of -^30 discrep- 
ancy. Indeed, the fraction of "missing" absorbers required would 



be inconsistent with established limits on the number of obscured 
DLA systems at high redshift (Ellison et al. 2001, 2004). 

Although the contribution of the low- and high-redshift Ca 11 
samples to the global SFR density is constant, the significant differ- 
ence between the observed SFR per unit area for the two samples 
is something of a puzzle. A resolution to the difference awaits the 
results of imaging studies of the Ca 11 absorbers and we will report 
on such observations in the near future. 

The low measured SFR associated with the Mg Il-selected 
DLA absorbers results in new constraints on the SFR per unit area 
and, assuming the absorbers are large (> 9h^^ kpc), on the contri- 
bution to the volume averaged SFR density, p* , of the host galaxy 
light falling within a ~ 7.5h~^ kpc radius of the DLA absorbers. If 
the absorbers are indeed large, and assuming the Rao et al. (2006) 
mean value for A^(H I), then the SFR per unit area may be as much 
as an order of magnitude, and certainly a factor of several, below 
that predicted by the Schmidt law. Correspondingly, the contribu- 
tion ofDLAs to p* (at 2~1) of the universe is dramatically smaller 
than that deduced by Hopkins et al. (2005), who assumed that the 
Schmidt law does apply. 

It is possible to make our constraints on the SFR per Mg 11- 
selected DLA candidate consistent with the predictions of the 
Schmidt law, and hence produce a contribution to p* as derived by 
Hopkins et al. (2005), but only by making the individual DLA ab- 
sorbers very small. Specifically, to bring the SFR per vmit area for 
the Mg ii-selected DLA candidate absorbers into agreement with 
the Schmidt prediction and {N(R I)) of Rao et al. (2006), the phys- 
ical size of the absorbers must be ~ 3(4)/t^^ kpc in radius, where 
the absorber size comes about by requiring that the surface area of 
an absorber times the SFR per unit area predicted by the Schmidt 
law, equals the integrated SFR per absorber we observe (Table 1). 
The size of the actual regions associated with star formation (i.e. 
with iV(H I)> 5 X 10^" cm~\ is then only ~ 2(2.5)/i~^ kpc in 
radius. The small physical size of the absorbers then requires that 
the space density of DLA absorbers at z^l is very large indeed, 
scenario in which such small absorbers are 
hosted by galaxies of similar space density is inconsistent with the 
established observation that the bulk of star formation at 2~1 oc- 
curs in luminous galaxies (Hippelein et al. 2003; Wolf et al. 2005; 
Bell et al. 2005). Alternatively, gaseous galaxy halos may contain 
many small DLA absorbers (i.e. have "filling factors" considerably 
less than unity), extending far beyond the extent of the SDSS fi- 
bre. Such a scenario is, however, inconsistent with observed DLA 
and associated metal line profiles as many independent components 
would be visible within the profiles (Wolfe, Gawiser, & Prochaska 
2005). In summary, reducing the size of the absorbers such that 
Schmidt law holds, produces a size/space-density constraint incon- 
sistent with any current favoured model for DLA absorbers and 
galaxies in general. 

Alternatively, the low-redshift Schmidt law may not hold for 
gas contained in intermediate to high redshift quasar absorption line 
systems. Recently, Wolfe & Chen (2006) have deduced an upper 
limit to the SFR per unit area associated with DLA absorbers at 
redshifts z = 2.5 — 3. 5 of less than 10% of that predicted from 
the Schmidt law. The Wolfe & Chen limit is based on the non- 
detection of the absorbers in the Hubble Ultra-Deep Field and is 
therefore a strict upper limit. Wolfe & Chen's interpretation is that 
in galaxies with low dust and molecular content, as in the DLAs, 
the threshold H I column density required to trigger star formation 
is significantly higher than in the nearby galaxies where the empir- 
ical Schmidt law was calibrated. Both these results, and the results 
presented in this paper, are inconsistent with the SFR per unit area 
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derived from C II* by Wolfe et al. (2003); for a full discussion of 
possible reasons we refer the reader to Wolfe & Chen (2006). 

Whatever the physical explanation, it is now difficult to escape 
the conclusion that only a small fraction of the star formation rate 
seen directly in galaxy surveys at redshifts from z = 0.5 to 3.5 is 
apparently associated with the relatively large 6-9 h^^kpc region 
surrounding DLA absorbers. Evidently, the largest contribution to 
the DLA cross-section is from gas which is too diffuse to support 
high rates of star formation and metal production, thus explaining 
the generally low metallicities of most DLAs (e.g. Akerman et al. 
2005). However, there are two plausible models consistent with our 
results. Recent theoretical work by Johansson & Efstathiou (2006) 
favours extended distributions of DLA gas, with the cross section 
dominated by small galaxies, and their predicted mean SFR per 
DLA host is in very good agreement with the values measured here. 
Recent observational work at both low and high redshift support 
such a scenario: Chun et al. (2006) measure luminosities of low 
redshift DLA hosts of ^ O.IL*; and Kulkami et al. (2006) place 
limits on the SFR of high redshift DLA hosts of 0.9 - 2.7M0yr"\ 
Alternatively, Wolfe & Chen (2006) associate low surface density 
H I regions, in which a reduced star formation rate efficiency exists, 
with the periphery of luminous, star-forming galaxies atz ~ 3. The 
low redshift Schmidt law may indeed still apply in the inner regions 
of these galaxies but the typical host to DLA absorber separation 
must then exceed the ~ 7.5h~^ kpc radius corresponding to the 
SDSS fibres. It should be relatively straightforward to assess the 
validity of these models at z < 1, where the putative star-forming 
galaxies associated, but not coincident, with the DLAs should be 
relatively easy to identify from ground-based images. We have be- 
gun such a study and will report on the results in the near future. 
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APPENDIX A: CORRECTING THE REDSHIFT PATH 
FOR THE EFFECT OF DUST ASSOCIATED WITH THE 

ABSORBERS 

The probability of finding an absorption line system intervening in 
a quasar spectrum depends on the total redshift path between the 
quasar and us, the signal-to-noise ratio (SNR) of the quasar spec- 
trum and the equivalent width of absorber that we are searching for. 
The effect of dust in an intervening absorption system is to reduce 
the SNR of the quasar spectrum, thus reducing the effective redshift 
path-length that each spectrum provides for finding an absorber of 
given equivalent width. The quasar may also be entirely removed 
from the sample as its magnitude falls below the survey flux limit. 

To estimate the total effect of a population of absorbers with 
known distributions of line equivalent width, dust content and red- 
shift, on the background quasar population of our survey, we sim- 
ply imagine placing each absorber in front of every one of the input 
quasars in our survey. We calculate the redshift path available to 
find a dust-free absorber: 

{Az)j,no-dust = ^ X P{Zabs,j,Wj,mk) (Al) 

fc,mf.<19.1 

and that available to find a dusty absorber causing A magnitudes of 
extinction in the observed frame i band: 

{Az)j^dust = ^ AXfc X P{zabs:j, Wj,mk + Aj) 

k,mk+Aj<19.l 

(A2) 

where the sum is over all quasars brighter than our survey mag- 
nitude limit, k denotes a quantity associated with the background 
quasar and j a quantity associated with the absorber to be found. 



AX is the redshift path defined by the Umits of the survey given in 
Section 2.1: 

AX = min[1.3, z^so - 0.1] - maa;[0.4, ^^^"^TOe"'"'"'' " ^' 

(A3) 

The probability, P, of detection of a line of given equivalent 
width and redshift, in a spectrum of a quasar of given magnitude, 
is calculated from Monte Carlo simulations: Gaussian line features 
are placed in the quasar spectra, our initial detection algorithm is 
run and the rate of recovery of the lines is measured. Including the 
probability of detection of the Mg II AA2796,2803 doublets has a 
negligible effect on the final results, due to their being considerably 
stronger than the Ca ii lines (mean equivalent width of Wa2796 = 
2.16A). 

The extinction in individual absorbers is estimated in the same 
way as for the subsamples, except that the absorber spectra are 
not combined into composites before the extinction curve is fit- 
ted. While the error on each individual reddening measure is large 
due to the unknown colour of the background quasar, the observed 
trends with equivalent width of Ca II give us confidence that on 
average the correct value is measured. 

A correction factor, Azno-dust/ Azduat, is calculated for 
each observed absorption line system, which tells us the number 
of absorbers missed for each one seen. An estimate is thus ob- 
tained of the total number of absorbers missing from our sample. 
Of course, this method only allows us to estimate correction fac- 
tors for regions of absorber parameter space which are well sam- 
pled by our observations. At high reddening values the number of 
absorbers drops dramatically and we are unable to correct reliably 
for absorbers with an E(B—V)> 0.2 and > 0.25 for the low- and 
high-redshift samples respectively. Our results are therefore only 
applicable to absorbers with reddening values below these limits. 
To place these E(B—V) limits into context, the reddening present 
in galaxies at z ~ 2, selected via ultraviolet continuum flux, has a 
median of 0.15 mag and ranges between 0.0 and 0.4 mag (Shapley 
et al. 2005). 
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